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ABSTRACT 


This dissertation presents a static analysis of the problem of an 
elastic layer perfectly bonded, except for a frictionless interface 
crack, to a dissimilar elastic half-plane. The free surface of the 
layer is loaded by a finite pressure distribution directly over the 
crack. The problem is formulated using the two dimensional linear 
elasticity equations. Using Fourier transforms, the governing equations 
are converted to a pair of coupled singular integral equations. The 
integral equations are reduced to a set of simultaneous algebraic 
equations by expanding the unknown functions in a series of Jacobi 
polynomials and then evaluating the singular Cauchy— type integrals. 

The resulting equations are found to be ill-conditioned and, 
consequently, are solved in the least— squares sense. 

Results from the analysis show that, under a normal pressure 
distribution on the free surface of the layer and depending on the 
combination of geometric and material parameters, the ends of the 
crack can open. The resulting stresses at the crack-tips are singular, 
implying that crack growth is possible. The extent of the opening and 
the crack-tip stress intensity factors depend on the width of the 
pressure distribution zone, the layer thickness, and the relative 
material properties of the layer and half-plane. 


v 




CHAPTER I 


INTRODUCTION 

Because of their high specific modulus and strength, advanced 
composite material systems have the potential to reduce the weight of 
aircraft structures. However, as with any new material system, the 
mechanical behavior of the material must be understood before it can be 
used extensively in structures. Understanding the mechanical behavior 
of composite materials is very challenging because of the complexity of 
the interactions between fiber and matrix, and between individual plies 
in a multilayered configuration. With composites have come not only the 
opportunity for the designer to "tailor" the material to optimize the 
structure, but also the challenge of a set of potential problems which 
were largely unknown in metal structures. For example, most of the 
composite systems are brittle, have low failing strains compared with 
metals, are susceptible to foreign object impact damage, and can develop 
delaminations (separation between plies). 

Some aspects of the delamination problem are addressed in the pres- 
ent work. Delaminations can be caused by manufacturing deficiencies, 
standard service loads, or extrinsic loads such as foreign object 
impacts. The damage which develops from impacts is a complex network of 
cracked plies, delamination between plies, and broken fibers; and it is 
very difficult to detect even with ultrasonic or radiographic techniques. 
A thorough understanding of impact damage is essential, as even small 
amounts of damage can substantially reduce both the tensile and compres- 
sive strength of a composite structure (see, for example, Rhodes [1]). 



Considering the delamination problem from the viewpoint of struc- 
tural life, it must be determined whether the delaminations are likely 
to grow under subsequent loading and may therefore place limits on 
either the loads or the life. The problem to be analyzed herein is 
chosen because it represents a basic mode of failure and should provide 
an indication of whether a delamination in layered composite materials 
might grow under the influence of purely compressive loading normal to 
the delamination. It is felt that a complete solution to the problem 
using the two-dimensional linear elasticity equations can contribute to 
a fundamental understanding of the mechanics of delaminations in layered 
media. Of course, seeking a complete solution necessitates restrictions 
on the complexity of the geometry and the constitutive equations. In 
contrast, more complex geometries and loadings can be analyzed using 
approximate methods such as finite element analysis. However, the power 
of the stress singularity at the ends of the delamination, a quantity 
critical to the application of elastic fracture mechanics methods, and 
other fundamental information on the behavior of the structure would be 
difficult or impossible to extract from a finite element analysis. The 
complete solution gives the form of the stress singularity directly from 
the governing equations. The present solution involves both the con- 
cepts of elastic contact problems and the analysis techniques associated 
with the singular stresses at the ends of the delamination (the delami- 
nation is modeled as a crack). 

Elasticity solutions of many contact problems and fracture problems 
involving flaws or cracks have been reported in the literature. For a 
comprehensive review of contact problems, the reader is referred to 
Gladwell’s recent book [2]. Of particular importance to the present 



investigation are papers by Erdogan and his colleagues [3, 4, 5, 6] and 
Keer et al. [7], in which the frictionless contact between elastic 
layers and elastic or rigid foundations is studied. All of the above 
solutions are for either two-dimensional or three-dimensional axisym- 
raetric problems. The layer is of finite thickness and the elastic 
foundation is either a half-plane or a half-space. The layer is taken 
to be weightless in all of these studies except in [4] and [5]. The 
layer rests on the foundation and load is applied to the free surface of 
the layer either as a compressive normal distributed load or through a 
stamp. In all of these problems, near the loaded region there is con- 
tact between the layer and the foundation. At some distance from the 
» 

loaded region (on the order of the width of the load or stamp) the layer 
separates from the foundation and comes back into contact with the 
foundation only when the weight of the layer is taken into account [5]. 
The contact stress is typically found to reach a maximum near the center 
of the loading and to vanish at the ends of the contact zone. The peak 
contact stress, the width of the contact zone, and the contact stress 
distribution depend on the layer thickness and the material properties 
of the layer and foundation. A somewhat different contact problem is 
solved by Keer and Chantaramungkorn [8]. In this problem, an elastic 
layer resting on an elastic half-space is loaded by a uniform normal 
compressive stress over the entire length of the layer except for a 
finite strip. The important conclusion of the study is that the layer 
separates from the half-space under the unloaded part of the layer. 

In all of these problems there is no bond between the layer and 
foundation. Thus, the question arises, "How will bonding or partial 
bonding change the stresses on the interface?" As a prelude to the 



present study, the axisymmetric problem solved by Keer et al. [7] was 
re-examined, using a solution technique similar to that used in [6], to 
investigate the effect of a modified boundary condition. The physically 
artificial boundary condition of no separation over a specified radius 
with zero shear stress was imposed. The three-dimensional axisymmetric 
solution was then found for a normally loaded elastic layer in contact 
over a circular region of prescribed radius with an elastic half-space. 
The results, which are consistent with those in [7], showed that for an 
arbitrary radius the contact stress was, in general, singular at the 
edge of the contact zone. However, for certain values of the contact 
radius the stress was zero at the edge of the contact zone. For contact 
radii which were about equal to the pressure distribution radius, the 
contact stresses were compressive, and for greater contact radii, the 
contact stresses near the end of the contact zone were tensile. 

The presence of these regions of tensile stress in the modified 
problem of [7] and of separation zones [8] under compressive loads 
suggest that if a layer were subjected to a compressive normal load 
directly over an interface delamination between the layer and founda- 
tion, then the layer could separate from the foundation near the ends of 
the delamination. Of course, crack opening might be expected only for 
certain combinations of layer thickness, delamination length, load 
distribution width, and material properties. It is known that for an 
opened crack, the crack-tip stresses are tensile and singular and hence 
the possibility of crack growth, leading to further weakening of the 
structure, exists. Typically, the magnitude of the stress intensity 
factors is related to the likelihood of crack growth. 



In [9] and [10], the behavior of layered elastic media with flaws, 
which are idealized as cracks, is studied. A uniform normal compressive 
stress is applied to both faces of the crack so that the crack is always 
open. Erdogan [10] finds the well known results that when the materials 
in the layers are identical the stresses exhibit a square-root singu- 
larity, and when the materials are different the stresses exhibit an 
oscillating square-root singularity. Stress intensity factors are found 
to depend on layer thickness and material property ratios. One should 
note that the oscillating singularity, which changes sign an infinite 
number of times within a small neighborhood of the ends of the crack, is 
physically unrealistic; yet it is the singularity dictated by the mathe- 
matics. The oscillating part of the singularity can be eliminated by 
modifying the boundary conditions near the ends of the crack [11]. The 
present problem, however, is formulated with the oscillating singu- 
^• ar liy> as this does not strongly influence the results away from the 
crack-tips. The feature of primary interest in the present contact 
problem is the separation point. Typically, the separation point is 
outside the region influenced by the oscillatory behavior. 

The problem analyzed herein is basically the same problem as that 
studied in [10] except for the external loading. The geometry is shown 
in Figure (la). However, changing the loading complicates the solution 
considerably. The problem is that of an elastic layer perfectly bonded, 
except for a crack on the interface, to a dissimilar elastic half- 
plane. Directly over the crack, uniform normal pressure is applied to 
the free surface of the layer. The problem is formulated as a two- 
dimensional elasticity problem; the solution is obtained using the 
singular integral equation techniques of Erdogan [12] . The integral 




a. Partial-contact. 



Figure 1. Geometry of the contact problems. 



equations derived herein are of the same fundamental form as those in 
[10] J the difficulty caused by changing the loading shows up in the more 
complicated boundary conditions in the contact region and in the fact 
that the extent of the contact region is unknown. The effects of these 
additional complications on obtaining a solution are discussed later. 

The solution gives contact stresses and displacements on the crack 
faces as well as crack-tip stress intensity factors. The results are 
presented as functions of material property ratios, the ratio of layer 
thickness to crack length, and the ratio of pressure distribution width 
to crack length. 



CHAPTER II 


FORMULATION 

The problem described in the Introduction is formulated using the 
two-dimensional linear elasticity equations. The solution is obtained 
by the use of Fourier transforms and singular integral equation tech- 
niques as described in [12] . The geometry of the problem is shown in 
Figure (la). The layer of thickness h and the half-plane, both isotro- 
pic, are assumed to be perfectly bonded along the interface except for a 
crack of length 2a. The contact between the crack faces is assumed to 
be frictionless. The load applied to the free surface of the layer is a 
uniform normal compressive stress, P Q , of width 2c located symmetrically 
with respect to the crack.. The region of primary interest in this 
problem is the interface. The solution gives stresses and displacements 
on the interface as well as crack-tip stress intensity factors. 

The governing equilibrium equations for the layer (i=l) and the 
half plane (i=2) are 


3a 


XX 


i 


+ 



and 


( 1 ) 


3a . 




3a 


yy, 


3x 


3y 


= 0 


( 2 ) 


where a xx ^ and Oyy are the normal stresses in the x and y directions, 
respectively, and a xy ^ is the shear stress. Using the plane strain 
stress-displacement relations 



9 


3u, 9v 

a xx i + 2u P ~bT + X i ly" * 


(3) 


3v. 3u. 

a yy ± " (x i + 2u i ) 1T + x ia3T ’ and 


/^i + M 

M i\ax 3y /» 


(4) 


(5) 


the equilibrium equations become 


w i v?u i + a i + ^ ib( 
“i^i + a i + ^ ir(' 


3u i 

— + — =°, and 


du 3v \ 

+ -=-i| = 0 


3x 3y 


( 6 ) 


(7) 


where X^ and are the Lame constants. Multiplying Equation (6) by 
sin(ax) and Equation (7) by 1/2/ ir cos( ox) and integrating both from 
zero to infinity yields the transformed equilibrium equations 


d 2 u. 


~(X i + 2p i )a 2 u i + — — - (X ± + u i )a = 0 , and 

dy ^ 


d2 ^i 2- du ^ 

(X ± + 2^) 2~ “ v i a v i + (x i + »*i>a 0 

dy 7 


( 8 ) 


(9) 


where ^ and v i are the Fourier transforms of u i and v ± defined by 

j~2 f- 

«i(o»y) ■ u ^( x »y) sin ox dx , and (10) 


-Jf f 


cos ax dx . 


( 11 ) 
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Simultaneous solution of Equations (8) and (9) yields 


= A ie “ oy + B^e ^ + Cf"’ + D^e^ , and 


-ay 


ay 




( 12 ) 


'i - v' ay + »i(<i W)*' w - v“ y + v“ y fi i 


(13) 


+ 3p i 


5A i + 6y i 


where — — for plane strain or for plane 

stress. The eight arbitrary constants (A^, Bj., C^, and D.^) are obtained 
from the boundary conditions which for the present problem are: 


cr (x,h) = 0 , 


(14) 


a (x,h) = P 

yy. 


,00 .["■ |x|<c , 

5 lo, |x| > c 


(15) 


a xy 1 ( x »°) = ^(x.O) , 


(16) 


a yyi (x,0) = a yy2 (x,0) ’ 


(17) 


u l( x >0) = u 2 (x,0), |x| > a 

o X y x (x,0) = 0, |x| < a 


(18) 


v l(x,0) « v 2 (x,0), |x| < b and |x| > a] 

Oy yi (x,°) =0, b < | x | < a 


(19) 


u 2 (x,-») = 0 , and 


( 20 ) 


v 2 (x,-») = 0 


( 21 ) 
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The contact distance b in Equation (19) is not known a priori but is 
found as part of the solution to the integral equations. The mixed 
boundary conditions, Equations (18) and (19), are difficult to utilize 
directly in developing a set of simultaneous equations for the eight 
arbitrary constants. Instead, it is convenient to define two unknown 
functions 

f l(x) = [u 1 (x,0 + ) - u 2 (x,0 )] , and (22) 


f 2 (x) = h [ v i( x * 0+ > “ v 2 ( x >° >1 


where the + and - superscripts refer to the limiting values of the dis- 
placements as y approaches zero from + and - sides, respectively • Using 
these two definitions in place of the mixed conditions, along with the 
remaining boundary conditions, Equations (14), (15), (16), (17), (20), 
and (21), the eight constants in the solutions for u^ and v^ are 
obtained* First, from Equations (20) and (21), which require that the 
displacements vanish for y + -<» y 

A 2 = B 2 = 0 ' 

The remaining six constants are obtained from the following six 
equations written in matrix form: 




where 
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^(cO 


f 2 (o) 





f^(x) cos(ax)dx 


f 2 (x) sin(ax)dx 


P Q (x) cos(ax)dx 


and 


(25) 


(26) 


(27) 


The solution to the six equations was obtained in closed form using 
MACSYMA, an algebraic manipulative computer code [13]. One should note 
that the mixed boundary conditions were not satisfied in obtaining the 
arbitrary constants. Instead, the constants are functions of f^ and f 2 , 
and the mixed boundary conditions are satisfied in formulating the 
governing integral equations. 

The solution to the problem is obtained by developing expressions 
for the stresses on the interface and then enforcing the mixed boundary 
conditions on those expressions. This yields two integral equations 
valid over (-a, a), with unknowns f^, f 2 , the contact stress (<Tyy(x,0), 

| x | < b) , and b. To obtain the integral equations for the interface 
stresses, only the solutions for C 2 and D 2 are required. They take 
the form 


C 2 oD(a) f[(c n ha+ c 12 )e a + (c 13 ha + c 14 )e ha ]P o e ha 
+ [c 21 + (c 22 h 2 a 2 + c 23 ha + c 24 )e 2ha + c 25 e -4ha ]f 2 

+ [ c 31 + ( c 32 h2a2 + c 33 ha + c 34 ) e 2h ° + c 35 e 4htX l^i ) * and 


(28) 
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D 2 - 5^7 ( [ <d U h “ + d 12 )e ' 2h '' + < d 13 h “ + d U> e ' 4h “l f o eh “ 

+ t d 21 + ( d 22 h2 “ 2 + d 23 ha + d 24 J e " 2h “ + d 25 e " 4h °! f 2 
+ [d 31 + (d 32 h 2 a 2 + d 33 ha + <* 34 )e _2ha + d 35 e" 4ha lf 1 } (29) 

where 

D(ot) = a^ + t a 12^ a ^ + a i3^ a + a iu) e + a 15 e . (30) 


and the c^j’s, dj-j's, and a i j f s are defined in Appendix A. In terms 
of u t and the transformed stresses are 

+ X 1 cS i , (31) 

(32) 

To obtain the expressions for the stresses on the interface, first 
evaluate a„„ and a T as y + 0 . Then, substituting u 9 and v 9 from 

yy2 x ^2 ^ ^ 

Equations (12) and (13) into Equations (31) and (32), using Equa- 
tions (28) and (29), and inverting (Xy , c xy 2 > anc * ^2 2 ives 


yyi 


V/ * J 0 yy ± 


cos( ax)dx = (X. + 2\i ) 


dv. 

l 


i i dy 


= J^I 0 Vj. sin<! “ )d>t ■ »i(dT - ®i)- 


Tia 


yy — 

= | J P Q Ct) ^ k^Cy.a) e ya cos( ax)cos( at) da dt 


[jfvjf 


+ / f~(t) I k 10 (y,a) e ya cos( ax)sin( at) da dt 

* / 0 L LZ (33) 


/"* i ( « /" 

•'O 1 * / 0 


k 13 (y,a) e^ a cos( ax)cos( at) da dt 


t _ , 

_ y=o 
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and 


7 TCT 

^2 y= 0 “ 

= I 

f" p # (t> I 

0 0 J 

r 09 

k 01 (y,a) e ya sin( ax)cos( at) 
0 Zi 

da dt 

S 

[J 



f°° 

r°° 



+ J 

f 2 (t) I 
0 J 

f k 22 (y»a) e ya sin( ax)sin(at) 

da dt 



f* 

C “ 



+ J 

f l(t) 

0 L J 

f k 9 o(y»a) e ya sin( ax)cos( at) 
0 

da dt 

where 





k n(y,a) = 

ha 

e 

D(a) 

{(l^jha + 

8 ^ 2 )e 2 * la + (g^ha + f^ 4 )e 



- y[d 

- 2 ha , 
12 ae + 

( d l 3 ha 2 + d 14 a)e 4ha ] } , 



(34) 


(35) 


- sby f B 21 + [ 8 22 h2 “ 2 + d 23 h “ + 8 24] e ’ 2h “ + S 25 e " 4h “ 


- y[ d 2l“ + ( d 2 2 h2 “ 3 + d 23 h “ Z + d 24“) e ' Z1, “ + d 25 ae_4ha ] ! . 


2 . , \_-2ha , , -4hai 


(36) 


k 13<y.«> - 5H0 (S 31 + («32 h2 “ 2 + 8 33 ha + 8,J<=" Zh “ + 8„»’ 4ha 


33 1 34 


'35 


- y[d 31 a + (d 32 h 2 « 3 + d^ha* + d 34 a)e Zha + d 35 ae~ 4ha ]} , 


-2ha 


-41101 


(37) 


ah 


k 2 i(y,a) - §^y {(0 41 ha + 6 42 )e 2ha + (B^ha + B 44 )e 


-4ha 


+ y[d 12 ae- 2hd + (d u ha 2 + d 14 «)e" 4h “] } , 


(38) 
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k 22<y>«> - D(a) ^51 + ( B 52 k a + %3 ka + %4 ) e + ^55 e 


'53 p 54 


55 


(39) 


+ y [ d 21 a + ( d 22 h2 ° 3 + d 23 ha2 + d 24 a ^ e ^ + d 25 ae > 


and 


k 23<y»«) D( a) ^61 + ^62 h a + B 63 ha + P 64^ e + S 65 e 


+ y[d-.a + (d 3 »h^a^ + d 3 ,ha^ + d»,a)« 


~2ha , , -4hcM i 
+ d 35 ae ]} . 


(40) 


The which are functions of material properties only, are defined 

in Appendix A. 

The fundamental form of the integral equations, Equations (33) and 
(34), depends on the nature of the integrands of the integrals with 
respect to a. Recognizing that some of the infinite integrals are not 
uniformly convergent for y = 0, the non-uniformly convergent parts of 
the integrands are separated from the convergent parts by writing 


k i 2 (y»«) = 


^ 21-^21 

a n a n 


ya - k’ 2 (y,a) , 


(41) 


k 13(y,a) = 


B 31 _ ^31 
a ll a ll 


ya - kj 3 (y,a) , 


(42) 


B 51 d 21 

k 2 2(y> a ) = * — + “ — ya - k' (y,a) , and 
a ll a ll Ll 


(43) 


k 2 3 (y,a) = ya - k' (y,a) 

a ll a ll 


( 44 ) 
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where the non-uni forraly convergent parts are the constant and linear 
terms in a and the uniformly convergent parts are the primed functions 
given by 



+ 





21 a. 


1_5 

11 > 


(45) 


4 3 (y,a) 





+ 


a. 


°31 a 


13 

11 



+ 



+ 



a 12V 2 2 

nrj h “ 


+ 



31 



+ 






» 


( 46 ) 
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The kernels k-^ and k£^ do not contain terras which are constant or lin- 
ear in a. The non-unif ormly convergent parts of the integrand, which 
can be evaluated in closed form, will provide the basis for determining 
the fundamental form of f^ and f£* Reference [14] evaluates infinite 
integrals such as those in Equations (33) and (34) with terras which are 
constant or linear in a. The results are reproduced in Appendix B. 
Using those results and Equations (41), (42), (43), and (44), the inte- 
gral equations. Equations (33) and (34), become 
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(X " p » (t) X -u 


* CO / k n (y>ct) e yct cos( ax)cos( at) da dt 
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C 00 
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-l 
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P 21 
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d 21 l \_ , + t 

3 11 [y 2 + (x + 


O 2 ] 2 [y 2 


x - t 

+ (x - t) 2 ] 2 


^2^’°^ e ya cos(ax)sin(at) da\dt 
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31 


2a 


11 


y 2 + (t - x) 2 


} 

+ ^— " J 

y + (t + x) z 


31 

2a y l 
11 


"y 2 ~ (t-x) 2 y 2 - ( t + x) 2 ~ 1 

_[y 2 + (t - x) 2 ] 2 [ y 2 + (t + x) 2 ] 2 J 


•'0 


k^Cy.a) e yct cos( ax)cos( at) da dt ^ dt] 


ly=0 


xy 2 _q- / / " / °° 

4^^ \J rt P o (t) J n e yot sin(ax)cos(at) da dt 


•jOK’ 


e, 


51 


y 2 + (t - x) 2 y 2 


1 

+ (t + x) 2 


21 


2a 


11 


2 

r V ~ 

_[ y2 + 


- (t - x) 2 _ y 2 - 

(t - X) 2 f [y 2 


■f; 


^22^*°^ sin( ax)sln( at) da 


- (t + x) 2 

+ (t +x) 2 ] 2 J 
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(49) 
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+ f f (t)/ ^ 61 f t + x t - x I 

•'O 1 V 2a H |_y 2 + (t + x) 2 y 2 + (t - x) 2 J 


+ -2i y 2 f-- + * t - x I 

a H [[y 2 + (t + X) 2 ] 2 [ y 2 +(t -x ) 2 ] 2j 

~ Jq ^23^ ,a ^ ® ya sin(ax)cos( at)daj dt^ J 


(50) 


The integrands of the infinite integrals with respect to a are of nega- 
tive exponential order for y < 0 and the limits exist. The limits of 
the remaining terms in Equations (49) and (50) are obtained from results 
in [14] (see Appendix B). Thus, the integral equations become 


ttct - 

yy 2 y=0~ f“ 

7"~ / P (t) I ^.(a) cos( ax)cos( at) da dt 

w 2 * / 0 ° * / 0 
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*2 W x + t x - 


dt 
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f 2 (t) / ^(a) cos( ax)sin( at)da dt 
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/ * 22 ^ sin (« x )sin(at)da dt 
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1-L J 


+ B -ii- « \ + a 15 

P 31 7^ 3 3A + e 3,, -X 3. 


31 a ll P35 / ’ 


cth _ 2 

*2l (ct) " nTTT t e a CBAi ha + 3 /l0 ) + 


-4ha, 


41 ha + 8 42 ) + e (3 43 ha + 8 44 >] , 



21 
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+ f 51 *n ' 654 ) 
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(57) 


i^3 (a) 


1 j -2ha 
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*u ' Ss3 ) 


ha 


+ f 61 a n " $64 J 


, -4ha/ 0 a 15 . 

6 (' 61 a ll " 65 J j 


(58) 


Noting from the boundary conditions that f^x) and f 2 (x) are zero for 
x •> a, that P Q (x) is zero for x > c and that f^(x) = f-^(-x), 

f2( x ^ = “ an( ^ P o< x > = p 0 (“ x )» t ^ ie integral equations can be 
written as 


■na 




y=0 


2 Vn 


f C f°° 671 f a f 2 (t > 

I P n ( t ) I <!-. (a) cos(ax)cos(at)da dt + I - dt 

J-c ° J 0 11 a ll J- a 1 ' x 


- /* h M f 

J-a J 0 

~P f 2 (C> / 

J -a * / 0 


63 1 

< 15 (a) cos(ax)sin(at)da dt + tt f. (x) 

a U 1 


^(t-) j <^3(0) cos(ox)cos(at)da dt, and (59) 


TTO 


yy 2 1 0 - fc s f" a., fa f,(t) 

^ “ | p (t) / KU-, (a) sin( ax)cos( at)da dt I — dt 

J-c ° J 0 a L1 J_ & t - x 


2y 


2 


-f-P^L 


Pc 1 

f 2 (t) ) <23 (a) sin(ox)cos(at)da dt + ir f 2 (x) 


'11 


- £ £ ^ (c) £ 


^(a) sin( ax)sin( at)da dt 


(60) 
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The equations for a yy^ an< * cr^^l q~ are va ^ ic ^ for an y x on y = 0 and 

the solution for f^(x) and f 2 (x) is obtained by restricting x to the 
interval (-a, a), i.e., the unbonded portion of the interface. 

In addition to the singular integral equations, the continuity 
conditions 


u i(x,0) = u 2 (x,0), |x| > a , and 

v l(x,0) = v 2 (x,0), | x | > a 


(61) 

(62) 


from the mixed boundary conditions. Equations (18) and (19), require 
that, in addition to f^(x) = = ® for l x l ^ a » 



dx = 0 , 


dx = 0 


and 


(63) 


(64) 


to prevent a rigid body displacement between the layer and the half- 
plane. Similarly, to prevent a relative rigid body displacement between 
the contact region |x| < b and the bonded region |x| > a, 


-b 

f 2 (x) dx = 0. (65) 

a 

Finally, the part of the mixed boundary condition, Equation (19), 
which requires continuity of normal displacement over the interior 
contact region, i.e., 



v l(x,0) = v 2 (x,0), 


I x | < b, 


( 66 ) 
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was not satisfied in formulating the integral equations and must be 
satisfied as an additional constraint on the solution to the integral 
equations . 

Equations (59) and (60), taken over (-a, a), and Equations (63), 
(64), (65), and (66) provide a complete set of equations to solve for 
the unknowns f^(x), f 2 (x), and °yy^ g* The remaining unknown, b, is 
obtained by an iterative process which is based on the requirement that 


the contact stress a. 


W2 


y=0 


vanish at x = ±b. The iterative process is 


described in Chapter III. 

Several fundamentally different problems can be solved using this 
set of equations. Two of these problems are described below. 

Full Contact in the Region (-a, a) ; 

For this problem, f 2 = 0 for all x and b = a, which places no 
restrictions on the sign of the contact stress. The integral equations 
reduce to a singular integral equation for f^ and a simple equation for 


yy 2 


y=0 


Equations (64) and (65) become trivial. The resulting equa- 


tions and the solution are presented in Appendix C. 


No Contact in the Region (-a, a): 

Two examples which involve no contact in the region (-a, a) arise 
(1) when the loading P Q (x) is prescribed as tensile, and (2) when the 
only prescribed loading is normal pressure on the crack faces, i.e., 


yy2 


is negative and <r is zero. For these problems the 

y=0 x ?2 y=0 


fundamental form of the integral equations does not change, but 

Equations (65) and (66) are no longer valid boundary conditions. The 

second problem (P-(x) = 0, a__, = 0, and a,.,. is compressive) 

° ^21 y=0 ^^2 1 y=0 

has been solved by Erdogan and Gupta [10] . 



CHAPTER III 


SOLUTION OF THE GOVERNING EQUATIONS 

Mixed boundary value problems such as the present 0 ne can usually be 
reduced to a system of singular integral equations with Cauchy-type 
kernels (such as Equations (59) and (60)). Erdogan et al. [12] give a 
detailed discussion of techniques for solving equations with 
Cauchy-type kernels. The present problem is solved using these 
techniques . 

Equations (59) and (60) are singular integral equations of the 
second kind as defined in [12]. By combining the equations into one 
complex integral equation, the singular behavior of fj. and f 2 can be 
determined from the dominant part of the integral equation. Then, 
expanding f^ and f 2 in a series of complex orthogonal polynomials, the 
singular terms in the integral equation can be removed and the integral 
equation can be reduced to a system of linear algebraic equations. The 
present problem has not only the unknowns f^ and f 2 , which appear in the 
singular integrals, but also o yy | y _ 0 and b. Further, the condition on 
normal displacement in the contact region. Equation (66), adds an addi- 
tional constraint on f^ and f 2 . The additional unknowns and the con- 
straint equation add considerable complexity to the solution of the 
problem. The details are discussed as the solution to the problem is 
developed. 

Defining r » ^ , s = — , and 
a a 


8l(s) = -fj/as) = — f L (x) , 


g 2 (s) = f 2 (as) = f 2 (x) 
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and requiring |x| < a, the Integral equations become 



K^Cct) cos(asa)cos(ta)da dt 


I P 8 2 <r) 
» J-l r - 


31 


B 


a, , 

, x , a 11 
8i(s) + “ R— 
1 tt 0 21 




k^Ccx) cos(asa)cos(ara)da dr 


a 11 


TT 0 


21 


IN (r) / 0 


cos(asa)sin(ara)da dr 


, U < Vy 2 << ‘ S 

2|1 2 B 21 




| s | <1 , and 


(67) 



i< 2 j(a) sin(asa)cos( ta)d 


1 p 8i<0 

a dt = — I - — 

* J-i r 


dr 


'51 , v . a a ll 


Hi g2<s> + * « 


61 


/ l 8 i (r> / 

J-l A 


< 23 ( 0 ) 


sin(asa)cos(ara)da dr 


a a ll 


ir B 


61 


/.N (r> Z" 


^ 2 ( 01 ) sin(asa)sin(ara)da dr, |s| < 1 


( 68 ) 


With the definition 


4>(s) = g 2 (s) + 1 gjjs) 


(69) 


the Integral equations combine to give 
- nr£ Ut - s)p o (t > 


H \ [Ki<r ' 


s)4>(r) + K 2 (r,s) <J>(r) ] dr 


+ i 6 2 a(s) 


|s| < 1 


y=0 


(70) 
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where 


L(t 


T h, 


(a) sin(aso)cos(ta) - i k^^(o) cos(asa)cos(ta) ] da , 


(71) 


K 


/ \ if"! /*22 (a) K i3( a ) 

l( r > 8 ' a j I ^ sin(asa)sin(ara) + — cos(asa)cos(ara) ] 

J 0 


) 


^l 2 (a) 
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+ il — -z cos(asa)sin(ara) — - — sin(asa)cos(ara)l 

\ P 21 %1 ) 


da, (72) 


K 2 (r,s) 


1 f°° r/ K i3 (a) 

T I cos(asa)cos(ara) - 

1_\ P 21 


<22^°^ 

B 61 


sin(asa)sin(ara)j 


^ 12 («) 


| iz 23 \ 

+ .11-7 cos(asa)sin(ara) H sin(asa)cos(ara)j 


<23(01) 


21 


61 


/ 


da , 
(73) 


<5i = 


B31 _ 8 51 ■ a n 


1 3 21 S 61 * ^ 2y 2 B 21 


and 


(74) 


CT < 8 > = ( as >l = ( as > I 


yy. 


2 I y=0 ^ 1 I y-0 


(75) 


Using Equation (69) the continuity Equations (63), (64), and (65), 
become 


/■ 


<f>(r) dr = 0 , and 


(76) 


ST 


Re(«(>(r)) dr = 0 , 


( 77 ) 


and the constraint equation. Equation (66), becomes 
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/ ReU(r)) dr = 0, - 7 < s < 7 . 

J- b/a a a 


(78) 


Because , K2, and L are sums of bounded functions, they are 
bounded as well. So, the singular behavior of <J> is determined from the 
dominant part of the integral equation, i.e., the first two terms on the 
right-hand side of Equation (70). From [12], this singular behavior has 
the form 


w(s) = (1 - s) 


“(1 + s) a , (w*(s) = (1 - s)“ a (l + s)" a ) 


(79) 


where 


1 — 1 1 1 ~ ^1 

a j “ i to, a = -j + ico, and co = ^ n \i~ T ~ 6" 


(80) 


Noting that w(s) is the weight of the Jacobi polynomials [12], it is 
natural to express the solution of the integral equation as 

CO 

<K S ) " w(s) C n P n a,a) ( s > • (81) 

n=0 


From the continuity condition. Equation (76), and the orthogonality 
relations [15] 

fo, n ^ m 



w(t) P^ ot » ct \t) P^ a,a \t) dx 
n m 


2 crt-od-1 

2n + a + a + 1 


x T(n + a + 1) T(n 4- a + 1) 
n! r(n + a + a + 1) 


( 82 ) 


n = ra , 


it is found that C Q = 0. 
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Substituting Equation (81) into Equation (70) and using the 
relation 




from [ 12 ], the integral equation reduces to 


(s) 


(83) 


-^jf 


L(t,s) P (t) dt 
c 0 n=l 


iU- 


- s 


- P ( “» a) ( 8 ) 


2 i x n -1 


V £1 [ K l (r ' s > C „ “ (r) P n°' 5> 


(r) 


+ K 2 (r,s) C n w(r) P 


n <a '“ ) ' r >] d r} 


+ i 62 a(s) , | s | < 1 


(84) 


where the singularity of the integral equation has been removed. 

Recognizing that for partial contact the normal stress on the crack 
must vanish at the ends of the contact region (s = ±b/a) the stress is 
represented as 



( 85 ) 


and the proper value of b/a is obtained by requiring that A Q = 0. 
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The unknown C n 's, A^'s, an< ^ b can be obtained numerically by 
reducing Equations (84), (77), and (78) to a set of linear algebraic 
equations. (The square set of equations generated initially from the 
governing equations was found to be ill-conditioned. Hence, to obtain a 
stable set of equations, the constraint equation was overspecified and 
the equations were solved in the least-squares sense.) 

The integral equation. Equation (84), can be reduced to a set of 
linear equations by expanding both sides of the equation as a series of 
Jacobi polynomials and then solving for the constant coefficients using 
the orthogonality relations. Equation (82). The resulting set of 
equations is 



Q (-a,-a) 

m 


C nri-1 


M 


+ E t B 


C + D cl 
n nm n J 


+ 


i 



R 

m 


> 


m =* 0,1 ,2,3, ... ,M . (86) 


where 


0 (-a,-a) 2 a rc*" 1 r(m - a + 1) r(m - q + 1) 

m 2m - a - a + 1 m! r(m - a - a + 1) * (87) 


a a 


11 


am 


/.I 


«*Cs) P<‘ a -’ 5) (S) 



P n“' S)<t) K l <r > s > 


dr ds 


> 


( 88 ) 
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D. 


nm 


it- x; "* (s) 
£ 5(r) 


(r) I^Crjs) dr ds , 


S £m 6 2 


A °f h/ «*( s ) P m" a> ° )(s) ds » 1 = 0 

l j ^ w*(s) P^" a ’“ a) (s) COB [|f2. (21 - lj ds , 


l - 1,2, • ,L , and 




T"Xi w * (s) P ® a ’~ a)(s) / c P o (t) L(t,8) 


dt ds . 


(89) 


(90) 


(91) 


The continuity equation. Equation (77) after substitution of 
Equation (81), can be integrated in closed form using the definition 
of the Jacobi polynomials [15] . The result is 


0 



r (- 1 ) 11 
n 2n 


(* - r 1 (> ♦ 


p (aU,etU)/b\ 
n-1 \a) 


(92) 


The constraint equation. Equation (78), also can be integrated in 
closed form and then both sides of the resulting equation can be 
expanded as series of Chebyshev polynomials [15] with constant coef- 
ficients. The coefficients are evaluated using the orthogonality 
relations for Chebyshev polynomials [15]. The resulting set of 
equations is 



( 93 ) 
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where 


- (1 - r) 0f< " 1 (1 + r )^ 1 pW.Wcrjj 


(94) 


and. rj is the Chebyshev polynomial of the second kind. 

The Equations (86) , (92) , and (93) yield M + Q + 1 linear algebraic 
equations for determination of M C n f s, and L + 1 A^'s. The C n f s and 
A^ f s are sufficient to determine the stresses and displacements along 
the entire interface, y = 0. Recognizing that the stresses on the 
interface are singular at s = ±1, the complex stress intensity factor 
is defined to be [10] 


k x + i k 2 


lim (x - 1) a (x + 1) a (a 
*♦!+ \ yy 2 


+ i 

y=0 



(95) 


Following [10] , the stress intensity factors are computed from 


ki + i k 2 


. + V j 1I 2 (I + 177 7 v- C 

C P|_ + l> 2 k 1^^2 + P 1 K 2^ V 1 2-# n n 

n=l 


(96) 


The set of simultaneous Equations (86), (92), and (93) must be 
solved numerically to obtain the C n 's and A^'s. As mentioned previ- 
ously, these equations were found to be ill-conditioned when L = Q [16] . 
Thus, a solution was sought by making Q > L and satisfying the equations 
in a least-squares sense [16]. Reasonable answers were then obtained. 
The condition of the set of equations was found to be further improved 
by setting M » L. This is consistent with the nature of the boundary 
conditions on the crack surface. The series on g 2 has M terras to 
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represent both the zero displacement difference in the region (-b/a,b/a) 
and the non-zero displacement difference in the regions (-1,-b/a) and 
(b/a,l); however, the stress, a(s) which has L terms, can be represented 
by only a few terms because it is smooth and symmetric. 

To obtain a valid solution to the problem, the value of b was 
incremented until |A q | < e where a satisfactory value of e depended on 
h. The solution was assumed to have converged when solutions with 
different values for M and L gave b f s which differed by less than 
one percent. 

The integrals in Equations (88), (89), (90), (91), and (93) were 
evaluated using quadrature formulas. Integrals with (1 - x 2)±l/2 
behavior near the end points were integrated using the Gauss-Chebyshev 
integration formulas [17]. The infinite integrals (K^(r,s) and K 2 (r,s)) 
were integrated using the Gauss-Legendre integration formula [17]. 



CHAPTER IV 


RESULTS AND DISCUSSION 
No Contact In the Region (-a, a) 

In this problem a crack of length 2a on the interface between the 
layer and the half-plane is opened by a unit normal pressure applied to 
both faces. The governing integral equations are given by Equa- 


tions (59) and (60) with P (x) = 0 and cr 


= -1. This problem, 


yy2| y =o 

which was first solved in [10], was re-solved herein to check the 
formulation and programming of the partial-contact problem. 

In [10] stress Intensity factors (for various values of h/2a) are 
presented for an epoxy layer and an aluminum half-plane (see Figure (6) 
in [10]). Corresponding results from the present formulation were found 
to be in excellent agreement with those results. 

A second check can be obtained by specifying identical elastic 
properties for the layer and the half-plane. For this problem the 
infinite integrals in Equations (59) and (60) can be evaluated in closed 
form. The resulting expressions are identical to those given by 
Equation (7.93) of [12], except for the sign of k 2 ^(x,t). It was 
verified in correspondence with Professor Erdogan that the minus sign 
had been omitted in [12] . 

The solutions of these two problems provide a check on the entire 
formulation of the partial-contact problem except for the terms 
involving P Q (x). They also validate the reduction of the integral 
equations to simultaneous algebraic equations as well as the programming 


of the solution. 
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Full Contact in the Region (-a, a) 

This problem contains the physically unrealistic boundary condition 
of frictionless adherence of the layer and the half-plane (see 
Figure (lb)). The resulting integral equations are considerably 
simplified, and they are readily solved numerically. Details of the 
formulation and solution of this problem are in Appendix C. This 
auxiliary problem was solved to obtain an initial estimate of the separ- 
ation point in the partial-contact problem. In fact, in all cases the 
point where the contact stress changed sign in the full-contact problem 
was an upper bound for the separation point in the partial— contact 
problem. 

Interestingly, the integral equations for the two problems differ 
considerably. The full-contact problem is governed by a single integral 
equation of the first kind, while the partial-contact problem is 
governed by a pair of integral equations of the second kind. In the 
full-contact problem, the stresses are singular at the crack-tip. The 
normal stress, a yy | y=0 , is singular as the crack-tip is approached from 
within the crack and non-singular as the crack-tip is approached from 
the bonded side. On the other hand, the shear stress, a^jy-Q, which 
is of course zero along the crack, is singular as the crack-tip is 
approached from the bonded side. This behavior of the stresses near a 
closed crack-tip has been discussed by Comninou [11]. In the partial- 
contact problem, the stresses have an oscillating square root singu- 
larity at the crack-tip. Both o yy (x,o) and o^yCxjO) are singular as the 
crack-tip is approached from the bonded side and, as the boundary 
conditions require, they are zero on the crack faces near the crack-tip. 
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The contact stresses on the crack are shown in Figure (2) for the 
case of a steel layer and an aluminum half-plane and in Figure (3) for 
the case of an aluminum layer and a steel half-plane. Both figures show 
normalized contact stresses versus position on the crack for several 
values of As is typical for contact between layers, the magnitude of 
the stress is largest for thin layers and it is concentrated under the 
applied pressure. For thicker layers the stress is distributed over a 
larger area. For the case of a steel layer and an aluminum half-plane 
(Figure (2)), the contact stress for the given geometries is compressive 
over the entire crack and at the crack-tip the stress is singular in 
compression. This behavior is typical for cases where > E2 and 

— > 0.3. When — and — are small the stress distribution can have a 

3 3 3 . 

region of tensile stress near the ends of the crack. However, the 
singular contact stress at the crack-tips is always compressive. For 
the case of an aluminum layer and a steel half-plane (Figure (3)), the 
contact stress is compressive over most of the crack and tensile near 
the crack-tip; the stress at the crack-tip is singular in tension. 

The point where the contact stress changes sign moves toward the crack- 
tip as the layer thickness increases. 

Figures (2) and (3) show that, if the modulus of the layer is 
greater than that of the half-plane, the singular stress is compressive, 
but that if this relationship is reversed, the singular stress is 
tensile. If the material properties are identical, the contact stress 
at the crack-tip is non-singular. This can be seen by examining the 
free term in Equation (C.l), which determines the contact-stress 
singularity. The material property coefficient, 6 31 , of the free term 
is zero when the material properties are identical. 
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Figure (4) shows the effect of the width c of the loaded region 

(relative to the crack length) on the contact stress; the contact stress 

is plotted versus position on the crack for several values of — . All 

a 

the results are for an aluminum layer and a steel half-plane. The 
magnitude of the contact stress Increases as the width of the loaded 
region increases, but the singular stress remains tensile, even for a 

c 

— value greater than one. 

. Figures (2), (3), and (4) establish several important trends in the 
behavior of the full-contact problem which can be carried forward to the 
discussion of the singular stresses at the crack-tips and later to the 
partial-contact problem: (1) the sign of the singular stress at the 

crack-tip depends on the relative stiffness of the layer and the half- 
plane; (2) changing the layer thickness to crack length ratio changes 
the contact stress distribution but not the sign of the singularity; and 
(3) changing the load width to crack length ratio changes the contact 
stress distribution but not the sign of the singularity. 

Using Equation (95), stress intensity factors for the material 
combinations of Figures (2) and (3) are plotted in Figures (5) and (6) 
as a function of dimensionless load width. One should note that k^ is 
not the classical mode I stress intensity factor but is the coefficient 
of the singular component of the contact stress on the unbonded side 
rather than the stress on the bonded side. The figure shows that for 
both material combinations and various values of — , the stress intensity 

cl 

£ 

factors increase to a maximum near — = 1.0 and then decay asymptotically 
to zero. Further, for both material combinations k£ is much greater in 
magnitude than k^. 
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The sign of reflects the sign of the singular stress at the 
crack-tips. For the aluminum layer and steel half-plane the singular 
stress is tensile and k^ is positive; for the steel layer and aluminum 
half-plane the stress is compressive and k^ is negative. This indicates 
that for a relatively "soft" (aluminum) layer the no-separation con- 
straint actually does prevent the layer from separating from the half- 
plane. For a relatively stiff layer the stresses at the crack-tips are 

He 

compressive, but for small — and — there can be a region of tensile 

3 . 3 

stress, which suggests the possibility of crack opening in the partial- 
contact problem. Interestingly, k 2 is negative for both material 
combinations, which indicates that the shear stress along the bond line 
constrains the layer from moving away from the origin relative to the 


half-plane, no matter which material is stiffer. 

max inside Vi 

With k^ and k 2 defined as the maximum k^ and k 2 for a given — 

(note that k^ x and k™ aX are functions of ■£•), Figure (7) shows k™ ax 

and k!? 331 plotted versus — for both material combinations. The peak 
ll 

values occur at — = 0*4 for the steel layer, aluminum half-plane 
a " 

h 

combination and at — 5 0.7 for the aluminum layer, steel half-plane 
combination. For large k™ ax and k^ aX become very small because of 
load diffusion effects. 


Figure (8) shows k^ 331 and k^ 3X plotted against the stiffness ratio, 
E 1 /E 2 for several values of X For E 1 /E 2 > 1, k“ ax is negative, and 
for E^/E 2 < 1, k^ ax is positive (Figure (8a)). This indicates that at 
the crack-tip the layer is attempting to pull away from the half-plane 
only when E 1 /E 2 < 1. The magnitude of kj 1331 is largest when the stiff- 


nesses are different, i.e., E^/E 2 - 7 or E^/E 2 - 0.1. The mode II 
stress intensity factor, k^ ax , is negative for all values of E 1 /E 2 . 



g- is not constant 


aluminum layer, steel half-plane 


steel layer, aluminum half-plane 


.2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 


Coefficient of the singular contact stress. 


is not constant 


/•^-aluminum layer, steel half-plane 
^ steel layer, aluminum half-plane 

. I I I I I I 1 1 I I 

.2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 


b. Mode II stress intensity factor. 


Figure 7. Variation of the maximum coefficients of the singular stress 
with layer thickness, — , for the full-contact problem. 
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For E 1 /E 2 + 0, k^ 3 * appears to approach a single value for all values of 
but that has not been established here. For E^/E 2 + <*>, kljf 3 * tends 
toward zero because the layer is rigid relative to the half-space. For 
an ideally rigid layer kj 13 * and k“ 3X are zero because a finite load on 
an infinite rigid layer will produce zero displacement at the interface. 

Partial Contact in the Region (-a t a) 

In this problem the layer is allowed to separate from the half- 
plane at the ends of the crack (see Figure (la)); the formulation and 
solution are in Chapters II and III. As discussed in Chapter II, the 
governing simultaneous equations are ill-conditioned and, consequently, 
solutions are obtained by satisfying the equations in the least-squares 
sense. The condition of the system of equations is found to depend on 
the combination of geometric parameters as well as the material 
properties of the layer and half-plane. As a result, satisfactory 
solutions are obtained for a relatively limited range of parameters, 
compared with the full-contact problem. To obtain solutions for a wider 
range of parameters would require more terms in the series expansion for 
4> (Equation (81)). However, the number of terms used to obtain the 
results presented herein required the maximum amount of storage 
available on the computer employed. 

Figure (9) shows the contact stress and normal displacement 
difference distributions for an aluminum layer and steel half-plane. 

The contact stress vanishes at about » 0.61. The normal displacement 
difference is essentially zero in the contact region, as required and 
the layer separates from the half-plane in the region of zero contact 
stress. The contact stress distribution is very similar to the 
compressive part of the distribution for the full-contact problem 
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(see Figure (3)). In the full-contact problem, the stress changes sign 
at about x ** 0.79, which is beyond the end of the contact zone in the 
partial-contact problem. The contact zone in the partial-contact 
problem should be smaller than the zone of compressive stress in the 
full-contact problem because, without the restraint imposed by the 
tensile stress of the full-contact problem, the crack will tend to open 
over an even greater area. Similar results were found for other 
combinations of geometric and material parameters. 

The dimensionless half-width of the contact zone, — , as well as the 

a 

magnitude of the contact stress and the crack opening, depend on the 

material and geometric parameters. Figure (10) shows the effect on — of 

a 

pressure distribution width, layer thickness, and material stiffness 
ratios. The size of the contact region increases with both layer 
thickness and pressure distribution width, but is insensitive to 
relative material stiffness for E^/E2 <1. As the ends of the contact 
zone near the crack-tips, substantial error develops In the displacement 
differences on the interface. For example, a negative normal displace- 
ment difference is predicted at the ends of the contact zone for an 
almost closed crack. Because the crack faces cannot overlap, this is 
physically impossible. The error develops for large contact regions 
because the large region of zero displacement difference prevents the 
series expansion for <|> from accurately modeling the small region of 
crack opening. For this reason, results are presented for relatively 

small — and — . 
a a 

The most important conclusion from the results for the partial- 
contact problem is that there is indeed crack opening near the ends of 
the crack under purely compressive external loading normal to the crack. 



a 


Effect of layer 
thickness. 


a 


b. Effect of pressure distri- 
bution zone width. 


£1 

e 2 

c. Effect of relative 
stiffness. 


Figure 10. Variation of contact zone width with geometric and material parameters. 


With crack opening comes the possibility of crack growth, and the 
likelihood of crack growth increases with the magnitude of the stress 

intensity factors. Figure (ll) gives kg as a function of h for several 

__ a 

values of c., and as a function of E^Eg. Figure (lla) shows that kg in 
a 

creases with both layer thickness and the width of the region of applied 
pressure. These results are consistent with those of the full-contact 
problem over the same range of parameters. Figure (lib) shows that kg 
increases with the material stiffness ratio E^/Eg. Again, these re- 
sults are consistent with results from the full-contact problem. Results 
for k^ are not presented because they appear to be of questionable 
accuracy. Computed values of k^ are typically an order of magnitude 
smaller than the values of kg, and thus are very sensitive to small 
errors in the series expansion of the displacement differences (which 
are used to compute k^ and kg). The apparent inaccuracy in k^ does not 
change the important result that the crack does open which is consistent 
with a positive k^. Because the normal displacement difference is 
closely related to k^ and it is constrained in the contact zone, k^ 
should be much smaller in magnitude than kg. This is in fact the case, 
although the values of k^ are suspect for the reasons mentioned above. 





51 


CHAPTER V 
CONCLUSIONS 


The problem of an elastic layer bonded, except for a frictionless 
crack, to a dissimilar elastic half-plane is considered in this study to 
investigate the behavior of the disbond under compressive loading. 

Using Fourier transforms, the two-dimensional elasticity equations are 
converted to a pair of coupled singular integral equations with Cauchy- 
type kernels. The equations are reduced to a set of simultaneous 
algebraic equations which are ill-conditioned; as a result, the equa- 
tions are solved in the least-squares sense. From the analysis, the 
following general conclusions about interface cracks in layered media 
are reached. 

1. Under a normal pressure distribution on the free surface of the 
layer the tips of the crack can open. The resulting stresses 
at the crack-tips are singular, implying that crack growth is 
possible even under compressive loading. 

2. The width of the crack-face contact zone depends on the layer 
thickness, the pressure distribution width, and the relative 
material properties of the layer and half-plane. 

3. Crack-tip stress intensity factors depend on the geometric 
and material parameters. 

4. Because the final set of simultaneous algebraic equations is 
ill-conditioned, accurate answers cannot be obtained for some 
values of the geometric and material parameters. To obtain 
accurate answers for a wider range of parameters, the program 
would have to be modified to allow for even more terms in the 
series expansion of the displacement-difference function. 
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Appendix A 

Material Property Terms 

The expressions for the material property terms in Equations (28), 
(29), and (30) are 

a ll " “(^<2 + ^2^**! + » 

a 12 = + ^2^ y l _ ^2^ * 

a 13 = 0 > 

a 14 " 2 * J i K 2 + y l y 2^ K l ” 1 ^ l< 2 - 1 ) “ V>2 ( K i + 1 ) » 
a 15 “ (^2^1 “ y l *2 ^ y l ” y 2^ * 

C 11 " \ + y 2^ K l + ^ » 

c 12 n 7 ” ^) 1 
c 13 " J C»*2 “ y l )K 2 (lc l + > 

c 14 = \ H 2 < k i + 1 )('< 1 - < 2 ) , 

C 21 " ” J ^(.< 1 ^ ~ » 


c 22 “ *" 2 y l ^2 ^ k 2 + l ) » 



c 23 " ^lH2 (lc l + 1)(k 2 “ 1 > » 
c 24 = J + x ) » 

c 25 " J \^ k 2 “ K i> * 

j- < 

C 31 * ^1^2 + V* [*1*2 + I (1 

c 32 = ^JV y 2^ K 2 ” ^ ~ ^ 1 * 2 ] 
c 33 = Wi^C^i + 1 )( |C 2 + 1) » 

2 1 

c 34 = ~ 2 »i *2 ~ 2 >*iH2^ ” 1)( 

c 35 3 h*2 " I v lV*l + k 2 ) ’ 

d ll a 0 » 

d 12 " ” J (^1 + li£ > C K x + 1) 
d 13 =* (U2 “ UiX^ + 1) * 

d 14 " 7 S " P2 )( *l + l) ’ 
d 21 “^(Ui + U2 k ].) » 

d 22 " ^(uj. ” \* 2 ^ » 


d 23 “ 2^ 1 ^( 1 ^ + 1) , 
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d 24 = 2 ^i[l J i + T " *)] » 

d 25 = C W2 “ V»i> » 
d 3i “ yj_ C p x + u 2 K i^ » 
d 32 ■ 4 Ui(p 2 ■ 1 * 1 ) * 

d 33 = 2 U 1 U 2 (k 1 + 1 ) , 

2 

d 34 “ “ 2 l J i “ “ 1 > » and 

d 35 ■ UiCWi - U 2 > * 

The definitions of the material property terms in the integral 
equations are 

e ll = -c ll » 

8 12 = J (1 + k 2> d i 2 - c u » 

813 " 2" (1 + k 2 ) d i3 - c i3 » 

®14 “ 2 ^ + k 2 * d 14 ** c 14 * 

®21 ” 2 ^ + d 21 ~ C 21 * 

e 22 = T (1 + k 2 ) d 22 " c 22 » 
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b 23 2 ^ + k 2 ^ d 23 “ c 23 * 

®24 7 ^ + k 2 ^ d 24 ” c 24 * 

025 ^ + k 2 ^ d 25 ~ c 25 * 

B31 ■ j (1 + k 2 ) d 3i “ c 3i » 

0 32 " J ^ + k 2 ^ d 32 ” c 32 * 

633 =■ ° , 

%4 “ 2 ^ + k 2^ d 34 ” c 34 * 

635 = ^ (1 + k 2 ) ~ , 

®41 C 11 ” 2 ^ k 2 ~ ^ d ll * 

®42 “ c 12 “ 2 ^*2 ~ ^ d 12 * 

643 — c^3 “ 2" ( k 2 — ^ d 13 » 

®44 a c 14 2 ^ k 2 ” d 14 * 

®51 “ C 21 ” 2 ^ k 2 ~ ^ d 21 * 

52 “ c 22 ” 2 ^ k 2 ” ^ d 22 * 
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6 53 = c 23 " \ (|C 2 " ^ d 23 » 

654 " C24 - ■g- (i<2 “ 1) » 

®55 ° c 25 ” 2 ^*2 ” ^ d 25 * 

%1 = C 31 “ 2 ^ K 2 “ ^ d 31 * 

®62 “ c 32 - 2" ^ K 2 ~ ^ d 32 * 

e 63 = c 33 “ J (k 2 " d 33 » 

6 64 “ c 34 ~ 2 ^ K 2 ~ ^ d 34 ’ 

%5 a c 35 ~ 2 ^ k 2 ~ ^ d 35 * 


and 
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Appendix B 

Relevant Integral and Limit Evaluations 


The following integral formulas, all for y < 0, are taken from [14]: 

I e yZ sin(x C )sin(tc)dc = j L XX - — = 1 T] 

J ° 2 [/ + (t - x ) 2 y 2 + (t + x) 2 J 


e y? x? sin(x?)sin(t?)d 5 = - X ■ y ~ ^ ^ 


2,2 


(y 4 + (t - x) 2 ) 2 (y 2 


- (t + x) 2 
+ (t + x) 2 ) 2 _ 


e a ^ sin(x£)cos(tc)dc = j 


1 t + x t - x 

2 Ly 2 + + x) 2 y 2 + (t - x) 2 _ 


/ e y? x ? sin(x ? )cos(tc)d ? = y 2 — * — — LJUE __T] 

* / ° (y + (t + x) 2 ) 2 (y 2 + (t - x) 2 ) 2 


cos(xc)cos(tc)dc - - 


e y ^ cos(x£)cos( t?)d £ = - 


2 [y 2 + (t - x) 2 y 2 + (t + x> 2 ] 


> and 


? 2 2 o “ 

X y - (t - x) y - (t + x) 

2 _(y 2 + (t - x) 2 ) 2 (y 2 + (t + x) 2 ) 2 


The following evaluations of limits of certain integrals are taken 


from [14]: 


_2 ± X. 

(t - x) 2 y 2 + (t 


dt = -m{»(y) , 
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Appendix C 

Solution of the Full Contact Problem 

For the full contact problem described in Chapter II, Equa- 
tions (59) and (60) reduce to 


ircx 


yy<j 


2p-j 


y=0 


£ £ 


V‘> J 0 ic^^(a) cos (at) cos (ax) da dt 


+ &31 f 

+ ir f 

a ll 


~ f^(t) J Kjj(a) cos Cot)ca8Cax)da dt , 


and 


(C.l) 


TTCT I 


CM 

If 

y=*0 


2 Ui 


£ p » (t> £ 
-/->>£ 


^(a) cos(at)sin(ax)da dt - 
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hi 


r a f x (t) 

J-a 


dt 


1^3 (a) cos(at)sin(ax)da dt 
The only continuity condition is Equation (63) , 

r a 

I f^(x) dx = 0 . 

J —a 


(C.2) 


(C.3) 


The only independent unknown in the problem is f ^ ; Oyy j is now 
a function of f^ . Thus, the solution for fj^ is obtained from Equa- 
tions (C.2) and (C.3). 

Restricting x to the interval (-a, a) and defining 


r “7 • 


s ■ — , and 
a 
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gj/s) - -f^as) - ~f j_(x) , 


Equations (C.2) and (C.3) become 


* P 61 J-c •'0 


P Q (t> f * 2 1 c°s(at)sin(asa)da dt 


i Pf iW 

* J-i r ~ 8 


dr 


a a l 1 f 1 /*" 

— J gj.(r) J ^(a) cos(ara)sin(asa)da dr , 


and 




g^(s) ds 3 0 . 


M < 1 . 

(C.4) 

(C.5) 


Equation (C.4) is a singular integral equation of the first kind 
(as defined in [12]), and the dominant part, 


0 61 g i (r) 


11 




dr 


determines the fundamental functional form of g^ which is 


gl(s) 


G 1 (s) 

Jl - s 2 


(C.6) 


where is a bounded continuous function. 

The numerical technique developed in [18] is used to solve Equa- 
tions (C.4) and (C*5). Following [18] the bounded function G^ is 
approximated as 


K 


G l( s) - £ C k I k (s) 
k-0 


(C.7) 
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where T k is the Chebyshev polynomial of the first kind. Then substi- 
tuting Equation (C.6) into Equation (C.4) with use being made of 
Equation (C.5), the integral equation reduces to the set of simultaneous 
equations, 


1 N x *1 

N £ W — - ~ V * k<s m- r n > * R(s m>- “ * 1 » " 

n-1 n m J 


(C.8) 


where use has been made of the appropriate relations from [18] and where 


a ll / 

k( r m » 8 n^ " a -»(<*) cos(ar a) sin(as a) da , 

P 61 *'0 


(C.9) 


1 a n f c r<* 

R ( s m ) “ " - J p c (t) cos(at)sin(as m a)da dt , 


n ■ cos ~ (2n - 1) , and 


(C.10) 

(C.ll) 


s 


m 


cos 



(C .12) 


Equation (C.8) gives N-1 equations for N unknowns; using the Gauss- 
Chebyshev integration formula [17], Equation (C.3) gives the N th 
equation as 

N 

W " 0 • (C.13) 

n»l 


After determining G x (r n ) by solving Equation (C.8) and C.13), the 
stresses on the entire interface are computed using Equations (C.l) and 
(C.2). The stress intensity factors are computed from Equation (95) 
with a = -1/2. 
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LIST OF SYMBOLS 


a half crack length 

a ij* c ij»^ij functions of material properties as defined In Appendix A 
b half width of the contact zone 


Ej.,E2 

f l. f 2 

81*82 

G 1 

h 

1 

kl»k2 

P 0 < x > 

P o 

p(a.a) 

n 

T k>Uk 

v l» u i 

w(x) 
a, a 

r 

6 1 ,6 2 

K i 


half width of the applied pressure region 

Young' 8 moduli of layer and half-plane, respectively 

displacement-difference functions on u and v, respectively 

dimensionless displacement-difference functions for f-^ and fo, 
respectively 

bounded continuous function defined in Equation (C.6) 
layer thickness 
/=! 

coefficients of the singular stresses 
pressure distribution on the free surface of the layer 
magnitude of the uniform applied pressure 
nth order Jacobi polynomial 

kth Chebyshev polynomial of the first and second kind, 
respectively 

displacements in the ith material in the y and x directions, 
respectively 

weight function of the Jacobi polynomials 
superscripts in the Jacobi polynomials 
material property terms defined in Appendix A 
gamma function 

material property coefficients as defined by Equations (74) 
and (75) 

material property constant for the ith material 



h>H 

o(s) 


Lame constants for the ith material 
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normal stress in the contact region 

'°w » a w stresses in the ith material 
**1 


"^2|jr-0 _ 

^ly-O- 

♦ 

0) 

V 2 


normal stress in the y-direction in the half-plane at y ■ 0 
shear stress in the half-plane at y = 0 
complex function of and g 2 defined by Equation (69) 
material property term defined by Equation (80) 

Laplace differential operator 


Subscripts 

x,y Cartesian coordinates 

1»2 layer and half-plane, respectively 
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